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Abstract. Based on the Nagel-Schreckenberg (NS) model with periodic boundary conditions, we proposed 
the NSOS model by adding the overtaking strategy (OS). In our model, overtaking vehicles are randomly 
selected with probability q at each time step, and the successful overtaking is determined by their velocities. 

We observed that (i) traffic jams still occur in the NSOS model; (ii) OS increases the traffic flow in the 
regime where the densities exceed the maximum flow density. We also studied the phase transition (from 
free flow phase to jammed phase) of the NSOS model by analyzing the overtaking success rate, order 
parameter, relaxation time and correlation function, respectively. It was shown that the NSOS model 
differs from the NS model mainly in the jammed regime, and the influence of OS on the transition density 
is dominated by the braking probability p. 

PACS. 89.40.-a Transportation - 45.70.Vn Granular models of complex systems; traffic flow 


1 Introduction 

Various dynamical models [1—4] have been proposed to 
explain the phenomena generated by traffic flow, such as 
fluid dynamical models, gas kinetic models and vehicle¬ 
following models [5-8] . From the microscopic point of view, 
vehicles can be represented by particles. Meanwhile, the 
way that vehicles influence other vehicles’ movements is 
treated as interactions among particles. The traffic system 
can be regarded to be composed of interacting particles 
far from equilibrium. Therefore, the vehicular traffic offers 
the possibility to study various fundamental dynamics of 
non-equilibrium systems, which are of interest in statisti¬ 
cal physics [2], 

During the last decades, cellular automata (CA) [9] 
have obtained popularity in simulating large networks, due 
to their simplicity and ability. One of the early CA based 
on traffic models is the NS model [ L0] developed by Nagel 
and Schreckenberg. Then, in order to make the CA more 
realistic, a large amount of modified models[ll-21] have 
been proposed by changing/adding some mechanisms on 
the NS model. Meanwhile, many theories and methods 
[22-30] that study the properties of the NS model have 
been proposed. For example, the jamming transition from 
a free flow phase to a congested phase is analyzed in de¬ 
tail by studying the relaxation into the steady state, mea¬ 
surements of an order parameter and the investigation of 
correlation functions [25-27]. 

The NS model is a probabilistic cellular automaton, 
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in which time and velocities are discrete, and the road 
is divided into sites. With periodic boundary conditions, 
the model deals with traffic flow of N vehicles moving in 
a one-dimensional lane of L sites. Each vehicle occupies 
one site. Each site can be either empty or occupied. Each 
vehicle has an integral velocity between 0 and v ma x (the 
maximum speed). Let v(j,t ) and x(J,t) denote the veloc¬ 
ity and position of the jth vehicle at time t, respectively. 
Let d(j, t) be the current empty sites in front of the j th 
vehicle, d(j,t) = x(j + 1 ,t) — x(j,t) — 1. The evolution 
rules of the NS model are as follows (parallel dynamics): 

(1) Acceleration, 

v{j,t i) -»• min(v(j,t) + 1 7 'Vmax)’i 

(2) Deterministic deceleration, 

v{j,t 2 ) ->■ min(v(j,ti),d(j,t)); 

(3) Random deceleration with probability p , 

v(j, t 3 ) -4- max(v(j , t 2 ) - 1,0); 

(4) Movement, 

x(j, t + 1) ->• x(j, t ) + v(j, t + 1). 

Here, ti, t 2l and t 3 are simply auxiliary intervals between t 
and t+1 for the convenience of implementing the CA rules 
computationally. The NS model is a very simple model 
which does not consider overtaking behavior. According 
to Step (2), the following vehicle is only able to move in 
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Fig. 1. (Color online) Illustration of updating and moving in 
the NS model and NSOS model, the red vehicle follows the 
blue one. (a) In the NS model, the red vehicle is only able to 
move in the regime d\. In the NSOS model, if the red vehicle 
is an overtaking one: (b) If the value of its velocity is smaller 
than d 2 , it follows the blue one. (c) If it could move more than 
dz, it overtakes the blue one. 
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Fig. 2. (Color online) Space-time diagram for L = 30, p = 0.3, 
q = 0.25, Vmax = 5 and p = 0.5. The Star stands for a free site, 
and the number stands for the velocity of a vehicle in that site 
after moving. With overtaking probability q, overtaking vehi¬ 
cles (colored rectangles) try to overtake the vehicle ahead, and 
the overtaking vehicles in the green rectangles have succeeded 
while the ones in the red rectangles have failed. Jams still occur 
in the case of p = 0.5. 


the d(j, 4) gap. For example, Fig. 1 (a) shows the mecha¬ 
nism of the original NS model, where the red car is only 
able to move within the gap d\ at time 4 + 1. However, 
in actual traffic situation, the red car could move much 
further within the gap d 2 (Fig. 1 (b)) or it could overtake 
the blue car, i.e., it could move further than the distance 
d 3 (Fig. 1 (c)). 

Considering the above facts, we proposed the NSOS 
model. In the real traffic, overtaking needs to be carried 
out on a surpass lane. If such a lane is occupied by other 
vehicles and not suitable for overtaking, the overtaking ve¬ 
hicle will become the ordinary one. Moreover, the human 
behavior is complex, and the overtaking decision is made 
by the driver. Due to these uncertainties, we introduce the 
probability q to describe the randomness of overtaking. In 
the NSOS model, each vehicle becomes an overtaking ve¬ 
hicle with probability q at each time step. Therefore, two 
kinds of vehicles exist in the NSOS model: one is the over¬ 
taking vehicle, and the other is the ordinary vehicle. Their 
types are varied according to the probability q. In fact, we 
ignore the detailed overtaking process and assume all the 
vehicles move on the same lane, based on the fact that the 
overtaking vehicles frequently return to the original lane. 

In this paper, we analyze the effects of overtaking by 
discussing the relaxation time, the spatial correlation func¬ 


tion and the order parameter. The results are compared 
to those of the NS model. Moreover, overtaking success 
rate we proposed was also studied, which can be used to 
depict the phase transition. 

The paper is organized as follows: In Section 2, the 
NSOS model is introduced. Simulation results are dis¬ 
cussed in Section 3, including the space-time diagram, fun¬ 
damental diagram, overtaking success rate, order param¬ 
eter, relaxation time and spatial correlation. Conclusions 
are given in Section 4. 


2 The NSOS model 

The NSOS model is a modified NS model incorporating 
overtaking strategy. In the NSOS model, every vehicle 
could be an overtaking one with probability q at each time 
step. This indicates that the system has qN overtaking ve¬ 
hicles on average at each time step. Overtaking vehicles 
will try to overtake the preceding ones, but whether they 
succeed or not depends on their own velocities and the 
positions of preceding vehicles at next time step. In or¬ 
der to avoid collisions (vehicle reaches the same position 
with the preceding one at next time step), the velocity 
of overtaking vehicle reduces by one. For simplicity, each 
overtaking vehicle is only able to overtake one vehicle per 
time step, and its position will locate in front of its preced¬ 
ing vehicle if it overtakes successfully. If two consecutive 
vehicles are both overtaking ones, and the preceding over¬ 
takes successfully, then the following one can not overtake 
successfully any more. The velocities of overtaking vehi¬ 
cles are reduced by one with probability p , except for the 
successfully overtaking vehicles. If the vehicle is an ordi¬ 
nary one, we update its velocity according to the rules 

(l)-(3) of NS model. The detailed updating rules are as 
follows: 

1 . At time 4, the jth vehicle becomes an overtaking vehicle 
with probability q, otherwise it is an ordinary one. 

2. Update the velocity: 

(I) If the jth vehicle is an ordinary one: 

(1) Acceleration: 

v(j,ti) ->• min(v(j, 4) + 1 

i V-max) • 

(2) Deceleration: 

v{j,t 2 ) -+ min{v(j,ti),d{j,t)). 

(3) Random braking: 

v(j,t 3 ) —>• max(v(j,t 2 ) — 1 , 0 ) with the proba¬ 
bility p. 

(II) If the jth vehicle is an overtaking one: 

(1) Acceleration: 

v(j,ti) -)• min(v(j, 4) + 1 

? V-max) • 

(2) If u(j, ti) > d(j , 4) + v(j + 1,4+1), the position 
d(j, 4) + v(j + 1,4 + 1) + 1 is empty and the 
(j + l)th vehicle does not overtake successfully, 

(i) Overtaking: 

v(j, < 3 ) -+ d(j, 4) + v(j + 1,4 + 1) + 1. 

(3) Otherwise, 

(i) Deceleration: 

v(j, t 2 ) -+ min(d(j, t)+v(j+l, t+l)-a, v(j, h)). 
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Fig. 3. (Color online) Fundamental diagrams for various values of overtaking probability q in the case of two random braking 
probabilities: (a) p = 0, (b) p = 0.25 with v ma x = 5. 
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Fig. 4. (Color online) Fundamental diagrams for various values of Vmax in different cases, (a) Deterministic NS model in the 
case of p = 0, q = 0. (b) NS model in the case of p = 0.25, q = 0. (c) NSOS model without random braking in the case of 
p = 0, q = 0.25. (d) NSOS model in the case of p = 0.25, q = 0.25. 


(ii) Random braking with probability p: 
v(j,t 3 ) ->■ min(v(j,t 2 ) - 1,0). 

3. Movement: 

x(j,t + 1) = x(j,t) +v(j,t 3 ). 

Here, v(j, t) denotes the velocity of the jth vehicle at time 
t and x(j,t) denotes its corresponding position. The num¬ 
ber of empty sites in front of the jth vehicle is denoted 
by d(j, t) = x(j + 1, t) — x(j, t) — 1. To avoid collisions, we 
assume a = 2 if the (j + l)th overtaking vehicle overtakes 
successfully. In other cases, a = 1. 

Since the velocity of overtaking vehicle at time t + 1 is 
relate to its preceding vehicle’s velocity, we should know 
the preceding vehicle’s velocity at time t + 1 first. For¬ 
tunately, the velocity of ordinary vehicle is independent 


of its preceding vehicle’s, so we could update the ordi¬ 
nary vehicles’ velocities first, and then update their rear 
overtaking vehicles’ velocities. In our simulations, we use 
the parallel update and periodic boundary conditions. In 
order to ensure the ordinary vehicles always exist in the 
NSOS model, which makes the update of overtaking vehi¬ 
cles work, we let the first and last vehicles be the ordinary 
vehicles all the time. 


3 Simulation results 

In this section, we will analyze simulation results of the 
NSOS model, including the space-time diagram, funda- 
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mental diagram, overtaking success rate, order parameter, 
relaxation time and spatial correlation. 

In our simulations, the system size is L = 1000 and the 
maximum velocity is assigned as v max = 5 unless noted 
otherwise. Simulations begin with vehicles randomly dis¬ 
tributed among the system. A single simulation runs for 
20000 time steps. The sampling for the calculation of the 
results is done for the last 10000 time steps and each data 
point is averaged over 100 different initial configurations. 

3.1 Space-time diagram 



P 

Fig. 5. (Color online) The overtaking success rate S' as a func¬ 
tion of density p for the three different conditions. The param¬ 
eters of the NSOS model are chosen as v max = 5 and p = 0.5. 



P 

Fig. 6. (Color online) The overtaking success rate S as a func¬ 
tion of density p for Vmax = 5, q = 0.5 and for different values 
of p. 

Fig. 2 presents a snapshot of a space-time diagram of 
the NSOS model simulated with L = 30, p = 0.3, q = 0.25, 
Vmax = 5 and p = 0.5. The space direction is horizontal, 
the time coordinate is downward, the number stands for 
the vehicle’s velocity after moving, and the star stands for 
a free site in which vehicle moves to the right. The colored 
rectangles contain overtaking vehicles. We can see that 
only a few (marked by green rectangles) overtake success¬ 
fully. In the red rectangles, unlike the ordinary vehicles, 
the overtaking vehicles could move further to the position 


that was occupied by their preceding ones at the previous 
time step, which is a common phenomenon on real traffic. 

We can see that congestion clusters still occur, and it 
indicates overtaking can not avoid traffic jams. The fact 
that sudden jams occur is not surprising since the NSOS 
model has much more situations which require strong brak¬ 
ing compared to the NS model. For example, in the top 
green rectangle in Fig. 2, the speed of overtaken vehicle 
changes from one to zero at the following step, which will 
induce a jam. 

3.2 Fundamental diagram 

In order to study the effects of the overtaking strategy in 
the NSOS model, we generate the fundamental diagrams 
under different traffic conditions. The average observable 
flow / versus density p is presented in the fundamental 
diagram. 

Firstly, we plot in Fig. 3 the fundamental diagrams 
with maximum velocity v max = 5 in various situations 
of different overtaking probability q and different random 
braking probability p. One observes that all the curves 
collapse in free-flow regime, since the overtaking vehicles 
make no contributions to the flow due to no interactions 
between them. It is interesting to notice that the values of 
maximum flow are suppressed for non-vanishing q. The 
major reason is that the overtaking mechanism makes 
more vehicles move closer, which will lead to more emer¬ 
gency brakes. Unlike the NS model, it is found that there 
exists a plateau at the large value of q. This is because 
the larger the overtaking probability, the more overtaking 
vehicles exist in the system that move synchronously with¬ 
out interacting with each other. Moreover, unlike in Fig. 3 
(b) where q = 0.25, the plateau never decays for p = 1 in 
Fig. 3 (a) where q = 0, as the absence of random braking. 
Another phenomenon is the flow enlarges as q increases in 
the jammed regime, which means that overtaking strategy 
increases the traffic flow in the jammed therein. 

Fig. 4 shows the fundamental diagrams for various val¬ 
ues of v m ax in different cases. Fig. 4 (a)-(d) are the deter¬ 
ministic NS model (p = 0, q = 0), the standard NS model 
(p = 0.25, q = 0), the deterministic NSOS model(p = 
0, q = 0.25), and the NSOS model ( p = 0.25, q = 0.25), re¬ 
spectively. Firstly we can see that the density correspond¬ 
ing to the maximum flow p(f ma x ) is apparently declined 
with the increase of v max ■ In fact, in comparison with 
Fig. 4 (a) and (b), asp increases, p(f max ) is also declined. 
While the value of p(f max ) does not change even when q is 
changed. Therefore, the main factors that determine the 
values of p{f max ) are the maximum velocity v max and the 
braking probability p, while the overtaking probability q 
just enlarges the flow in the jammed regime. 

Specifically, according to [25], only in the determinis¬ 
tic NS model (Fig. 4 (a)), the flow is solved exactly given 
by f{p) = min(pv max , 1 - p), and p(f max ) coincides with 
the critical density p c , which is the transition density from 
the free flow phase to the jammed phase. In other words, 
P^fmax) is not equal to p c except in the deterministic NS 
model. 
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3.3 Overtaking success rate 

Since not all the overtaking vehicles can overtake success¬ 
fully in the NSOS model, it is natural to explore the over¬ 
taking success rate. Here, an overtaking vehicle is marked 
by s, and s = 1 if it overtakes successfully, otherwise s = 0. 

The overtaking success rate per time step per vehicle 
is defined by 


t 0 +T Nq 


-| -i •'UT-c 

s = ?v; E £»*■ 


t=t 0 +1 2=1 


(1) 


where No is the number of overtaking vehicles in the sys¬ 
tem, to is the time sufficiently long after the system has 
reached a steady state, taken as to = 10 5 . 

In fact, overtaking success rate S' is a function of the 
density p, braking probability p and overtaking probabil¬ 
ity q. First, we show in Fig. 5 the overtaking success rate 
S as a function of density p calculated under different 
overtaking probabilities q. The simulation parameters are 
v max — 5 and p = 0.5. At low densities, overtaking will 
not occur until the density reaches a critical density p c . 
This is because in the free-flow region (p < p c ) all vehicles 
move with the velocity v max . The average gap between 
vehicles is superior to the maximum speed v max and thus 
successful overtaking vehicles do not exist. 

From Fig. 5, we observe that the values of the critical 
densities p c almost remain the same. This is because brak¬ 
ing probability p , which essentially determines the critical 
densities shown in Fig. 6, is the same in the three different 
conditions. 

Above the critical density, S increases with the in¬ 
crease of density, and reaches a maximum value, then de¬ 
creases with much larger density. Comparing the three 
different conditions, we can see that the larger overtaking 
probability leads to larger value of S. At relatively high 
density, the overtaking success rate S reaches the max¬ 
imum value. From our daily traffic experiences on high¬ 
ways, overtaking usually occurs when the vehicles are par¬ 
tially in jams and free flow. 

To show how S depends on the braking probability 
p , we plot in Fig. 6 the values of S against density p 
with maximum velocity v max = 5, overtaking probabil¬ 
ity q = 0.5 and for different values of p. We observe that 
the critical density p c decreases when increasing p. An¬ 
other interesting result is that beyond p c , the overtaking 
success rate S increases when p decreases. These results 
could be explained by the fact that when the braking prob¬ 
ability p increases, the vehicles’ movements become slow. 
In a certain way, the overtaking is more difficult. 


3.4 Order parameter 

In the NS model, the phase transition, from the free flow 
regime to the jammed phase, has been investigated [25- 
27]. Here, we use the same approaches to analyze the ef¬ 
fects of overtaking strategy in the NSOS model. 

To look for the transition between the free and jammed 



Fig. 7. (Color online) Order parameter m vs. density p (L = 
1000) for various values of q in the case of p = 0.25. They do 
not vanish exactly for p < p c but converges smoothly to zero. 



P 

Fig. 8. (Color online) The behaviors of order parameter m 
and overtaking success rate S near the critical density. Here, 
p = 0.25 and q = 0.25. 



P 

Fig. 9. (Color online) Relaxation time near the transition 
density for different system sizes for the deterministic case 
(p = 0,q = 0). 


phase, an order parameter has been introduced in [25], 
which is defined by 


1 ! t 0 +T L 

m= TL J2 riirli + 1 ’ ( 2 ) 

t=io + l i=1 

with rij = 0 if a site is empty, and n, = 1 for a site occu¬ 
pied by a vehicle. In fact, the order parameter m is count¬ 
ing occurrence of neighboring vehicles, which describes the 
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Fig. 10. (Color online) Relaxation time near the transition 
density in the NSOS model and NS model under p = 0.05. 

jamming level. 

In the determined NS model {p = 0), the transition 
density is p c = l/(v m ax + 1)- Below the transition density, 
the values of order parameter are zero because every vehi¬ 
cle moves with v max . Within the jammed phase, the flow 
is limited by the number of empty sites and motionless 
vehicles occur. 

When p > 0, the behavior of the order parameter qual¬ 
itatively changes in the vicinity of the transition density 
in the NS model. It does not vanish exactly at p < p c but 
converges smoothly to zero even for small values of the 
braking probability p. 

In order to investigate the effects of overtaking in the 
NSOS model, we plot in Fig. 7 the values of to against the 
vehicle density p with maximum velocity v max = 5, brak¬ 
ing probability p = 0.25 and for various values of q. We 
find that the order parameters in the NSOS model have 
the similar behavior in the NS model, which belongs to 
the continuous phase transition. The disparity is that the 
values of order parameters are slightly increased with the 
growth of the values of q in the jammed phase. This result 
could be explained by the fact that the NSOS model pro¬ 
vides much more configurations of adjacent vehicles than 
the NS model, which leads to the increase of order param¬ 
eter to. Again, we also observe that the critical density 
remains unchanged with respect to the different overtak¬ 
ing probabilities. 

Fig. 8 shows the behaviors of order parameter m and 
overtaking success rate S near the critical density. One 
can see that around the critical density the order param¬ 
eters are the same as the overtaking success rate S, which 
further confirms that the overtaking success rate is able 
to identify the critical state. In fact, below the critical 
density, there is no difference between the order parame¬ 
ter and overtaking success rate, while the disparity occurs 
above the critical density. 


3.5 Relaxation time 

A characteristic feature of a second order phase transition 
is the divergence of the relaxation time at the transition 
point. According to Ref. [26], the relaxation time is charac¬ 
terized by the parameter (Due to the dimensionless units 



r 

Fig. 11. (Color online) (a) Correlation function in the vicinity 
of the phase transition for the deterministic case. At p = p c 
the amplitude is independent of the distance r. (b) In the case 
of q = 0, the amplitude of the correlation function decays ex¬ 
ponentially. (c) In the case of q = 0.1, the amplitude of the 
correlation function decays exponentially faster compared with 
(b). 

no normalization is introduced) 

[min{v*(t), (Foo)} - {v(t))]dt, (3) 

where v(t) is the average velocity of all vehicles as a func¬ 
tion of time t starting from randomly positioned vehicles 
with zero initial velocity (v(t = 0) = 0). For t —» oo the 
system reaches a stationary state with average velocity 
(voo). v*[t) denotes the average velocity in the acceleration 
phase t —>• 0 for low vehicle density p —>• 0, v*(t) = (1 —p)t 
holds in this regime since the vehicles don’t interact with 
each other. Thus the relaxation time is obtained by sum¬ 
ming up the deviations of the average velocity (v(t)) from 
the values of a system with one single vehicle that can 
move without interactions with other vehicles (p —» 0). 
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In our simulations, the total running time of system is 
10L + T. We calculate the velocities when t < 10L. After 
10-L time steps, the system reaches the stationary state, 
and then we average the stationary velocity Vqo during 
the T time steps, with T = 5000 here. Each data point is 
averaged over 100 different initial configurations. 

At first, we have reviewed the scaling behavior of the 
relaxation time in the deterministic model shown in Fig. 9. 
We consider different system sizes up to L = 20000 where 
the position of the critical density p c becomes size inde¬ 
pendent. It is found that there exists a maximum value of 
relaxation time at the critical density. Taking the magni¬ 
tude of r as a characteristic value for the relaxation time, 
we can estimate the dynamical exponent. Moreover the 
scaling behavior of the width <j(L) and height r m (L) of 
the peak has the following form, 

T m (L) oc L z ,a(L) oc (4) 

We obtain the critical exponents for the deterministic case 
(p = 0, q = 0) are z = 0.53 and p = 6.7. 

In order to study the effect of overtaking on the relax¬ 
ation time r, we assign p = 0.05 and q = 0.1 in the NSOS 
model, and make a comparison with the NS model (q = 0). 
The results are shown in Fig. 10. Unlike the deterministic 
case, the shape is not symmetric. We find that the relax¬ 
ation time of the NSOS model (q = 0.1) is larger than 
that of the NS model (q = 0) at critical density, which 
implies that the system will take more time to reach the 
stable state due to the extra fluctuations induced by the 
overtaking probability q. Another phenomenon is that the 
negative relaxation times occur both in the two models, 
this is because the system gets temporarily into states that 
have a higher average velocity than the stationary state, 
i.e., (v(t)) > (Too), in the case of p > p c [25]. During the 
first few time steps, small clusters occur in the initial con¬ 
figuration, and then more and more vehicles get trapped 
into large jams. Therefore, the average flow decreases to 
its stationary value. 


3.6 Spatial correlation 

Another feature of second order phase transitions is the 
occurrence of a divergence length scale at critical den¬ 
sity and a corresponding algebraic decay of the correlation 
function. The correlation function is defined as 



Fig. 12. (Color online) Correlation length according to the 
maximum amplitude of correlation function in different q in 
the case of p = 0.008 at critical density p c = 0.1666. The dots 
represent the results of simulations, and the solid lines are ex¬ 
ponential fittings. The fitting parameters are the corresponding 
correlation length. For the NS model (q = 0), its correlation 
length is £ = 25.5. For the NSOS model in the case of q = 0.1, 
q = 0.3 and q = 0.5, respectively, their correlation lengths are 
£ = 10.8, 4.5 and 3.5. 
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Fig. 13. (Color online) Density dependence of the correlation 
length in the vicinity of the phase transition density with p = 
0.008, q = 0. 
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g{r) = —j ’Y^mni+r - p 2 , (5) Fig. 14. (Color online) T'h.e ninximum correlation length ^max 

^ L t=t 0 +1 2 —i as a function of the overtaking probability q for p = 0.008. 


where n t = 0 if the lattice is empty, otherwise n, = 1. 
Here, we obtain the same result in the vicinity of the phase 
transition for the deterministic NS model (p = 0, q = 0) 
shown in Fig. 11 (a). There is a decay of the amplitude of 
g(r) for larger value of the distance between the sites. Pre¬ 
cisely at pc, the correlation function is solved exactly given 
by g{r) = p c - p 2 c for r = 0 rnod(v max + 1) or g(r) = -p 2 c 
for otherwise. 


In order to study the effect of overtaking on the cor¬ 
relation function, we have compared the results of the 
NS model with NSOS model at small value of p, here 
p = 0.008, shown in Fig. 11 (b) and Fig. 11 (c). The com¬ 
mon feature is that their amplitudes of the correlation 
function decay exponentially, except that in the NSOS 
model such a decay is much faster. The reason is that in 
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the NS model after the transient state, the vehicles arrange 
themselves in regular intervals and therefore the ordering 
is long-range. In the NSOS model, with random braking 
and/or random overtaking, this long-range ordering is de¬ 
stroyed. 

Furthermore, we calculate the correlation lengths in 
the two models according to their maximum amplitude 
shown in Fig. 12. For the NS model (q = 0), the corre¬ 
lation length is £ = 25.5. While in the NSOS model, the 
correlation length is shorter, for example, £ equals 10.8 
when q = 0.1. These results demonstrate that the overtak¬ 
ing reduces the spatial correlation. The maximum value of 
correlation length £ max determines the critical density for 
small value of p in the NS model. For example, Fig. 13 
shows the relationship between correlation against den¬ 
sity in the case of p = 0.008 and q = 0. We could observe 
that the correlation length reaches the maximum value 
at the critical density. We also calculate £max for differ¬ 
ent overtaking probabilities q shown in Fig. 14. We can 
see that £ max does not always decay with the growth of 
overtaking probability q , and there is a critical value of 
q corresponding to the minimum correlation length. Here, 
in the case of p = 0.008, the critical overtaking probability 
is about q c = 0.5, and the minimum correlation length is 
£max — 3.5. 


4 Conclusions 

This paper has proposed the NSOS model for a simplified 
traffic situation. Based on the NS model, we add the mech¬ 
anism of overtaking strategy by bringing in the overtaking 
probability q. Basically, one major difference between the 
two models is that our model makes the overtaking pos¬ 
sible, allowing one vehicle overtakes the preceding one or 
moves further within larger gap. 

We have studied the properties of NSOS model un¬ 
der periodic boundary conditions, and made a compari¬ 
son with the results of NS model. The space-time diagram 
tells us that the traffic jams are still spontaneous forma¬ 
tion, even if the existence of overtaking vehicles. And we 
find traffic flow enlarged in the jammed regime, according 
to the comparison of the fundamental diagrams with that 
of NS model. Furthermore, we redo many phase transition 
calculations within the NS framework, such as the order 
parameter, relaxation time and spatial correlation. It is 
found that the critical density is dominated by the braking 
probability p other than overtaking probability q. Accord¬ 
ing to our model, we proposed the concept of overtaking 
success rate S. We find it could be used to characterize 
the phase transition, since successful overtaking vehicles 
do not occur in the free flow phase. 
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